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We present the first search for heavy, long-Uved particles that decay to photons at a hadron 
collider. We use a sample of 7+jet+missing transverse energy events in pp collisions at \/s = 
1.96 TeV taken with the CDF II detector. Candidate events are selected based on the arrival time 
of the photon at the detector. Using an integrated luminosity of 570 pb^^ of collision data, we 
observe 2 events, consistent with the background estimate of 1.3±0.7 events. While our search 
strategy does not rely on model-specific dynamics, we set cross section limits in a supersymmetric 
model with Xi and place the world-best 95% CL. lower limit on the x? mass of 101 GeV/c^ 

at r^o = 5 ns. 

PACS numbers: 13.85.Rm, 12.60.Jv, 13.85.Qk, 14.80.Ly 
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Searches for events with final state photons and miss- 
ing transverse energy {frr) [H at coUider experiments 
are sensitive to new physics from a wide variety of mod- 
els [2^] including gauge mediated supersymmetry breaking 
(GMSB) In these models the lightest neutrahno (x?) 
decays into a photon (7) and a weakly interacting, stable 
gravitino (G) that gives rise to fjr by leaving the detec- 
tor without depositing any energy. The observation of an 
ee^^fjT candidate event by the CDF experiment during 
Run I at the Fermilab Tevatron [3] has increased the in- 
terest in experimental tests of this class of theories. Most 
subsequent searches have focused on promptly produced 
photons 0, 0| , however the Xi can have a lifetime on the 
order of nanoseconds or more. This is the first search 
for heavy, long-lived particles that decay to photons at a 
hadron collider. 

We optimize our selection requirements using a GMSB 
model with a standard choice of parameters [7| and vary 
the values of the x? mass and lifetime. However, the final 
search strategy is chosen to be sufficiently general and 
independent of the specific GMSB model dynamics to 
yield results that are approximately valid for any model 
producing the same reconstructed final state topology 
and kinematics § . In pp collisions at the Tevatron the 
inclusive GMSB production cross section is dominated 
by pair production of gauginos. The gauginos decay 
promptly, resulting in a pair of long-lived Xi's in asso- 
ciation with other final state particles that can be identi- 
fied as jets. For a heavy x? decaying inside the detector, 
the photon can arrive at the face of the detector with a 
time delay relative to promptly produced photons. To 
have good sensitivity for nanosecond-lifetime x^'s jS], we 
search for events that contain a time-delayed photon, 
and > 1 jet. This is equivalent to requiring that at least 
one of the long-lived Xi's decays inside the detector. 

This Letter summarizes the first search for heavy, 
long-lived particles that decay to photons at a hadron 
collider. The data comprise 570±34pb~^ oipp collisions 
collected with the CDF II detector [iBl at Vs = 1.96 TcV. 
Previous searches for nanosecond-lifetime particles using 
non-timing techniques yielded null results 11 1 . 

A full description of the CDF II detector can be found 
elsewhere ^ . Here we briefiy describe the aspects of the 
detector relevant to this analysis. The magnetic spec- 
trometer consists of tracking devices inside a 3-m diame- 
ter, 5-m long superconducting solenoid magnet that op- 
erates at 1.4 T. An eight-layer silicon microstrip detector 
array and a 3.1-m long drift chamber with 96 layers of 
sense wires measure the position {xi) and time [ti) of the 
pp interaction [l^ and the momenta of charged particles. 
Muons from collisions or cosmic rays are identified by a 



system of drift chambers situated outside the calorime- 
ters in the region with pseudorapidity |?7| < 1.1 [ij. 
The calorimeter consists of projective towers with elec- 
tromagnetic and hadronic compartments. It is divided 
into a central barrel that surrounds the solenoid coil 
(1 77 1 < 1.1) and a pair of end-plugs that cover the region 
1.1 < \r]\ < 3.6. Both calorimeters are used to identify 
and measure the energy and position of photons, elec- 
trons, jets, and fr- The electromagnetic calorimeters 
were recently instrumented with a new system, the EM- 
Timing system (completed in Fall 2004) , that mea- 
sures the arrival time of electrons and photons in each 
tower with \r]\ < 2.1 for all energies above ~5 GeV. 

The time and position of arrival of the photon at the 
calorimeter, tf and Xf, are used to separate the photons 
from the decays of heavy, long-lived Xi'^ from promptly 
produced photons or photons from non-collision sources. 
We define the corrected arrival time of the photon as 



^University of California Irvine, ''IFIC(CSIC-Universitat de Valen- 
cia), 



t2 = tf-- u 



The tj distribution for promptly produced, high energy 
photons is Gaussian with a mean of zero by construction 
and with a standard deviation that depends only on the 
measurement resolution assuming that the pp production 
vertex has been correctly identified. Photons from heavy, 
long-lived particles can have arrival times that are many 
standard deviations larger than zero. 

The analysis preselection is summarized in Table HI It 
begins with events passing an online, three-level trigger 
by having a photon candidate in the region \r]\ < 1.1 
with Et> 25 GeV and |t> 25 GeV. Offline, the high- 
est Et photon candidate in the fiducial region of the 
calorimeter is required to have Et > 30 GeV and to 
pass the standard photon identification requirements [H| 
with a minor modification [3]. We require the event 
to have > 30 GeV where the trigger is 100% effi- 
cient. We require at least one jet with |?7j°*| < 2.0 and 
S^"' >30 GeVlH. Since a second photon can be identi- 
fied as a jet, the analysis is sensitive to signatures where 
one or both Xi's decay inside the detector. To ensure 
a high quality ti and Xi measurement, we require a ver- 
tex with at least 4 tracks, X^tracks^'T > GeV/c, and 
\zi\ < 60 cm; this also helps to reduce non-collision back- 
grounds. For events with multiple reconstructed vertices, 
we pick the vertex with the highest X^tracksP^- 
duce cosmic ray background, events are rejected if there 
are hits in a muon chamber that are not matched to any 
track and are within 30° of the photon. After the above 
requirements there are 11,932 events in the data sample. 

There are two major classes of background events: col- 
lision and non-collision photon candidates. Collision pho- 
tons are presumed to come from standard model interac- 
tions, e.g., 7-t-jet-(-mismeasured fir, dijet-|-mismeasured 
^ where the jet is mis-identified as a 7, and W ^ ev 
where the electron is mis-identified as a 7. Non-collision 



5 



Preselection Requirements Cumulative (individual) 

Efficiency (%) 

> 30 GeV, It > 30 GeV 54 (54) 
Photon ID and fiducial, jjyj < 1.0 39 (74)* 
Good vertex, EtracksP^ > 15 GeV/c 31 (79) 
|7?j<='| < 2.0, Sj."* > 30 GeV 24 (77) 
Cosmic ray rejection 23 (98)* 
Requirements after Optimization 
^ > 40 GeV, > 35 GeV 21 (92) 
A(?!)(|?r, jet) > 1 rad 18 (86) 
2 ns < t? < 10 ns 6 (33) 

TABLE I: The data selection criteria and the cumulative 
and individual requirement efficiencies for an example GMSB 
model point at m^o = 100 GeV/c^ and r^o = 5 ns. The ef- 
ficiencies listed are, in general, model-dependent and have a 
fractional uncertainty of 10%. Model-independent efficiencies 
are indicated with an asterisk. The collision fiducial require- 
ment of Izil < 60 cm is part of the good vertex requirement 
(95%) and is estimated from data. 



backgrounds come from cosmic rays and beam effects 
that can produce photon candidates, pr^ and sometimes 
the reconstructed jet. We separate data events as a func- 
tion of tj into several control regions that aUow us to 
estimate the number of background events in the final 
signal region by fitting to the data using collision and 
non-collision shape templates as shown in Fig. [TJ 

Collision photons are subdivided in__two subclasses: 
correct and incorrect vertex selection 



13| . An incorrect 



vertex can be selected when two or more collisions occur 
in one beam bunch crossing, making it possible that the 
highest reconstructed X^tracks Pt vertex does not produce 
the photon. While the fraction of events with incorrect 
vertices depends on the final event selection criteria, the 
tj distribution for each subclass is estimated separately 
using W ep data where the electron track is dropped 
from the vertexing. For events with a correctly associ- 
ated vertex, the t2 distribution is Gaussian and centered 
at zero with a standard deviation of 0.64 ns [l3'|. For 
those with an incorrectly selected vertex the t2 distribu- 
tion is also Gaussian with a standard deviation of 2.05 ns. 

The tj distributions for both non-collision backgrounds 
are estimated separately from data using events with 
no reconstructed tracks. Photon candidates from cos- 
mic rays are not correlated in time with collisions, and 
therefore their tj distribution is roughly flat. Beam halo 
photon candidates are produced by muons that origi- 
nate upstream of the detector (from the p direction) and 
travel through the calorimeter, typically depositing small 
amounts of energy. When the muon deposits significant 
energy in the EM calorimeter, it can be misidentified as a 
photon and cause These photons populate predomi- 
nantly the negative region, but can contribute to the 
signal region. Since beam halo muons travel parallel to 
the beam line, these events can be separated from cosmic 



ray events by identifying the small energy deposited in 
the calorimeter towers along the beam halo muon trajec- 
tory. 

The background prediction uses control regions out- 
side the signal time window but well within the 132 ns 
time window that the calorimeter uses to measure the 
energy. The non-collision background templates are nor- 
malized to match the number of events in two time win- 
dows: a beam halo-dominated window at {—20, —6} ns, 
selected to be 3cr away from the wrong vertex collision 
background, and a cosmic rays-dominated window at 
{25, 90} ns, well away from the standard model and 
beam halo contributions. The collision background is 
estimated by fitting events in the {—10, 1.2} ns window 
with the non-collision contribution subtracted and with 
the fraction of correct to incorrect vertex events allowed 
to vary. In this way the background for the signal region 
is entirely estimated from data samples. The systematic 
uncertainty on the background estimate is dominated by 
our ability to calibrate the mean of the t2 distribution 
for prompt photons. We find a variation of 200 ps on 
the mean and 20 ps on the standard deviation of the dis- 
tribution by considering various possible event selection 
criteria. These contribute to the systematic uncertainty 
of the collision background estimate in the signal region 
and are added in quadrature with the statistical uncer- 
tainties of the final fit procedure. 

We estimate the sensitivity to heavy, long-lived parti- 
cles that decay to photons using GMSB models for dif- 
ferent Xi masses and lifetimes. Events from all SUSY 
processes are simulated with the pythia Monte Carlo 
program [3] along with the detector simulation [l3| . The 
acceptance is the ratio of simulated events that pass all 
the requirements to all events produced. It is used in 
the optimization procedure and in the final limit setting 
and depends on a number of effects. The fraction of x? 
decays in the detector volume is the dominant effect on 
the acceptance. For a given lifetime this depends on the 
boost of the Xi- A highly boosted Xi that decays in 
the detector typically does not contribute to the accep- 
tance because it tends to produce a photon traveling in 
the same direction as the Xi- Thus, the photon's arrival 
time is indistinguishable from promptly produced pho- 
tons. At small boosts the decay is more likely to happen 
inside the detector, and the decay angle is more likely 
to be large, which translates into a larger delay for the 
photon. The fraction of events with a delayed photon ar- 
rival time initially rises as a function of Xi lifetime, but 
falls as the fraction of Xi's decaying outside the detector 
begins to dominates. In the Xi mass region considered 
(65 < m^o < 150 GeV/c^), the acceptance peaks at a 
lifetime of around 5 ns. The acceptance also depends on 
the mass as the boost effects are mitigated by the ability 
to produce high energy photons or ^ in the collision, as 
discussed in Ref. Q. 

The total systematic uncertainty of 10% on the ac- 
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Standard Model 
Beam Halo 
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-20 20 40 60 80 
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FIG. 1: The time distribution for photons passing all but 
the final timing requirement for the background predic- 
tions, data, and a GMSB signal for an example point at 
m^o = 100 GeV/c^, r^o = 5 ns. A total of 1.3±0.7 back- 
ground events are predicted and 2 (marked with a star) are 
observed in the signal region of 2 < < 10 ns. 



ceptance is dominated by the uncertainty on the mean 
of the t2 distribution (7%) and on the photon ID effi- 
ciency (5%). Other significant contributions come from 
uncertainties on initial and final state radiation (3%), jet 
energy measurement (3%), and the parton distribution 
functions (1%). 

We determine the kinematic and t2 selection require- 
ments that define the final data sample by optimizing 
the expected cross section limit without looking at the 
data in the signal region. To compute the expected 95% 
confidence level (C.L.) cross section upper limit we 
combine the predicted GMSB signal and background esti- 
mates with the systematic uncertainties using a Bayesian 
method with a flat prior [19]. The expected limits are op- 
timized by simultaneously varying the selection require- 
ments for pr, photon Et, jet Et, azimuth angle be- 
tween the leading jet and frr {Aip{frr, jet)), and i^- The 
A(j}{f}r, jet) requirement rejects events where the fJr is 
overestimated because of a poorly measured jet. While 
each point in Xi hfetime vs. mass space gives a slightly 
different optimization, we choose a single set of require- 
ments because it simplifies the final analysis, while only 
causing a small loss of sensitivity. The optimized require- 
ments are summarized in Table HI As an example, the ac- 
ceptance for m-^o = 100 GeV/c^ and lifetime t^o = 5 ns 
is estimated to be (6.3±0.6)%. 

After all kinematic requirements, 508 events are ob- 
served in the data before the final signal region time re- 
quirement. Their time distribution is shown in Fig. [T] 
Our fit to the data outside the signal region predicts total 
backgrounds of 6.2±3.5 from cosmic rays, 6.8±4.9 from 
beam halo background sources, and the rest from the 




6 70 75 80 85 90 95 100 105 110 
mass (GeV/c^) 

FIG. 2: The contours of constant 95% C.L. upper cross section 
limits for a GMSB model 0. 



standard model. Inside the signal time region, {2, 10} ns, 
we predict 1.25±0.66 events: 0.71±0.60 from standard 
model, 0.46±0.26 from cosmic rays, and 0.07±0.05 from 
beam halo. Two events are observed in the data. Since 
the result is consistent with the no-signal hypothesis, we 
set limits on the Xi hfetime and mass. Figure[2]shows the 
contours of constant 95% C.L. cross section upper limit. 
Figure [3] shows the exclusion region at 95% C.L., along 
with the expected limit for comparison. This takes into 
account the predicted production cross section at next- 
to-leading order [2^ as well as the uncertainties on the 
parton distribution functions (6%) and the renormaliza- 
tion scale (2%). Since the number of observed events 
is above expectations, the observed limits are slightly 
worse than the expected limits. These limits extend at 
large masses beyond those of LEP searches using photon 
"pointing" methods |11|. 

In conclusion, we have performed the first search for 
heavy, long-lived particles that decay to photons at a 
hadron collider using data collected with the EMTim- 
ing system at the CDF II detector. There is no excess 
of events beyond expectations. As our search strategy 
does not rely on event properties specific solely to GMSB 
models, we can exclude any 7-l-jet-l-^ signal that would 
produce more than 5.5 events. We set cross section limits 
using a supersymmetric model with Xi ~^ jG, and find 
a GMSB exclusion region in the Xi hfetime vs. mass 
plane with the world-best 95% C.L. lower limit on the 
Xi mass of 101 GeV/c^ at t^o = 5 ns. Future improve- 
ments with similar techniques should also provide sen- 
sitivity to new particle decays with a delayed electron 
signature 0. By the end of Run II, an integrated lumi- 
nosity of 10 fb~^ is possible for which we estimate a mass 
reach of ~ 140 GeV/c^ at a lifetime of 5 ns. 
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Y+p^+1jet analysis with ElVITiming (570 pb'^) 
^^^1 Predicted exclusion region 
Observed exclusion region 
ALEPH exclusion upper limit 



GMSB x'^yG 
M,,,,=2A,tkn(P)=15 
N„„,,=1 , n>0 



is the azimuthal angle, and — — lntan(6'/2). The trans- 
verse energy and momentum are defined as Et = E sin 9 
and pt ~ psinO where E is the energy measured by 
the calorimeter and p the momentum measured in the 



tracking system, pr 
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~0 2 

mass (GeV/c ) 



FIG. 3: The exclusion region at 95% C.L. as a function of Xi 
lifetime and mass for a GMSB model The predicted and 
the observed regions are shown separately and are compared 
to the most stringent published limit from LEP searches [lH . 
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